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--■^ ■ 1. Introduction I: Remarks on chiral extrapolations 

|' The first part of this talk concerns the application of chiral perturbation 

theory (CHPT) at unphysical quark masses. More precisely, lattice QCD 
\ (LQCD) allows one in principle to calculate hadronic matrix elements ab 

initio using capability computing on a discretized space-time. To connect 
to the real world, various extrapolations are necessary: LQCD operates at a 
finite volume V, at a finite lattice spacing a and at large (unphysical) quark 
$J] ■ masses m q . All these effects can be treated in suitably tailored effective field 

theories (EFTs) (for a recent review, see 1 ). Here, I consider the quark mass 
expansion of certain baryon observables. Various nucleon (baryon) observ- 
ables have already been computed on the lattice, like e.g. masses of ground 
and excited states, magnetic moments, nucleon electromagnetic radii, the 
nucleon axial-vector coupling, and so on. CHPT in principle provides ex- 
trapolation functions for all these observables, parameterized in terms of 
a number of low-energy constants (LECs). These LECs relate many ob- 
servables, they are not dependent on the process one considers. Given the 
present situation with only a few lattice results at reasonably small quark 
masses available, it is mandatory to incorporate as much phenomenological 
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input as is available for the LECs from studies of pion-nucleon scattering, 
pion production and so on. Also, since many LECs appear in various observ- 
ables, a true check of our understanding of the chiral symmetry breaking 
of QCD requires global fits at sufficiently small quark masses. It is ab- 
solutely necessary for such extrapolations to make sense that one is in a 
regime where higher order terms stay sufficiently small. Consequently, re- 
sults must be independent of the regularization scheme, these can differ by 
higher order terms (e.g. comparing results based on heavy baryon CHPT to 
ones obtained employing e.g. infrared regularization). For this interplay of 
CHPT and LQCD to make sense, the lattice "data" should be in the true 
chiral regime. I will illustrate these issues for two very different examples: 
the axial- vector coupling gA and the Roper mass tur. 

2. Application I: The nucleon axial-vector coupling 

The nucleon axial-vector coupling gA is a fundamental quantity in hadron 
physics as it appears prominently in the Goldberger-Treiman relation. Lat- 
tice results for gA obtained by various collaborations for pion masses be- 
tween 300 and 1000 MeV show a very flat quark (pion) mass dependence. 
On the the other hand, it is known since long that the one-loop represen- 
tation of gA is not converging well and is dominated by the M 3 term with 
increasing pion mass, 2 thus gA{M n ) rises steeply as the pion mass increases. 
The large coefficient of this term is, however, understood in terms of the 
large values of the dimension two LECs C3 and C4 combined with some large 
numerical prefactors. In Ref. 3 we have therefore worked out the two-loop 
representation of gA, 



with g the chiral limit value of gA, A is the scale of dimensional regular- 
ization, and the coefficients 02,3, /?2 encode the one-loop result. Further, F 
denotes the chiral limit value of the pion decay constant, F ~ 86 MeV, and 
A(") collects the corrections ~ M™. At two-loop order, one has corrections 
of fourth and fifth order in the pion mass, given in terms of the coefficients 
0:4,5, (3a, 74 (note that there is also a M T In terms whose contribution 
we have absorbed in the uncertainty of 05). The LEC 04 can be analyzed 
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Fig. 1. The axial- vector coupling qa as a function of the pion mass. The solid (red), 
the dot-dashed (black) and the dashed (green) line correspond to various input values 
for the LECs (see 3 for details). The (magenta) circle denotes the physical value of g\ 
at the physical pion mass, the triangles are the lowest mass data from Ref. 8 and the 
inverted triangles are recent results from QCDSF. 9 

using the renormalization group (as stressed long ago by Weinberg 4 ) and 
is entirely given in terms of the dimension three coefficients of the one- 
loop generating functional. We have performed this calculation based on 
two existing versions of the dimension three pion-nucleon chiral Lagrangian 
without and with equation of motion terms (see 5 and, 6 respectively) and 
obtained 

16 11 o A 

a 4 = -y-y.g 2 + 16.g^ ■ (2) 

In 3 the dominant contributions to the LECs 0:5, ^4,74 from 1/rnjy correc- 
tions to dimension two and three insertions (with the nucleon mass) 
as well as from the pion mass expansion of F n were also worked out, for 
details see that paper. Setting the remaining contributions to zero, we find 
(in the notation of Eq. (1)) g A = g {l - 0.15 + 0.26 - 0.06 - 0.001) and 
go = 1.21, using gA = 1.267 and central values for the LECs C3, C4, die, see 
e.g. 7 This shows that for the physical pion mass, the higher order correc- 
tions are small and one thus has a convergent representation. Varying the 
LECs 0:5, /?4, 74 within bounds given by naturalness, one finds that the pion 
mass dependence of gA stays flat for M n < 350 MeV, see Fig. 1. From this 
figure one also sees that there is just a little overlap between the lattice re- 
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suits and the chiral representation, where it can be applied with a tolerable 
uncertainty. We remark that another solution to this problem was offered 
in Ref. 10 where an effective field theory with explicit delta degrees of free- 
dom at leading one-loop order could lead to a flat pion mass dependence of 
gA, requiring, however, a fine tuning of certain low-energy constants. For 
a recent update, see . n It should also be noted that most of the lattice 
results analyzed in these papers are far outside the range of applicability of 
that particular EFT evaluated only to leading one-loop order - that such 
a representation works at such large pion masses is an interesting obser- 
vation but certainly does not support claims of a controlled and precise 
determination of qa from LQCD. 

3. Application II: Chiral corrections to the Roper mass 

Understanding the (ir)regularitics of the light quark baryon spectrum poses 
an important challenge for lattice QCD. In particular, the first even-parity 
excited state of the nucleon, the Roper iV*(1440) (from here on called the 
Roper) is very intriguing — it is lighter than the first odd-parity nucleon 
excitation, the 5n(1535), and also has a significant branching ratio into 
two pions. Recent lattice studies have not offered a clear picture about the 
nucleon resonance spectrum. In particular, in Ref. 12 an indication of a rapid 
cross over of the first positive and negative excited nucleon states close to 
the chiral limit was reported - so far not seen in other simulations at higher 
quark masses. Note also that so far very simple chiral extrapolation func- 
tions have been employed in most approaches, e.g., a linear extrapolation 
in the quark masses, thus ~ M 2 , was applied in. 13 It is therefore impor- 
tant to provide the lattice practitioners with improved chiral extrapolation 
functions. A complete one-loop representation for the pion mass depen- 
dence of the Roper mass was recently given in. 14 Since the Roper is the 
first even-parity excited state of the nucleon, the construction of the chi- 
ral SU(2) effective Lagrangian follows standard procedures, see e.g. 15 The 
effective Lagrangian relevant for our calculation is (working in the isospin 
limit m u = and neglecting electromagnetism) 

C = Cq + Cr + Cnr , 
Co = iNj^D^N - m N NN + iR^D^R - m R RR , 

= \g R Ri^R , C% = i^R W u»N + h.c. , 
4 2) =cl{ x+ )RR- -^R^uJiD^D"} + h.c.) R+°-t(u^)RR , 

8777 n Z 
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= -^(x + ) 2 RR , (3) 

where N,R are nucleon and Roper fields, respectively, and toat,to_r the 
corresponding baryon masses in the chiral limit. The pion fields are col- 
lected in m m = —d^Tr/Fjr + 0(n 3 ). is the chiral covariant derivative, 
for our purpose we can set = cL, see e.g. 15 for definitions. Further, 
X+ is proportional to the pion mass and induces explicit chiral symme- 
try breaking, and ( } denotes the trace in flavor space. The dimension 
two and four LECs c* and e* correspond to the Cj and e, of the effective 
chiral pion-nucleon Lagrangian. The pion-Ropcr coupling is given to lead- 
ing chiral order by C R , with a coupling gR. This coupling is bounded by 
the nucleon axial coupling, \gR\ < \gA\, in what follows we use gR = 1. 
The leading interaction piece between nucleons and the Roper is given by 
C^' R . The coupling g^ r can be determined from the strong decays of the 
Roper resonance, its actual value is g^R — 0.35 using the Roper width 
extracted from the speed plot (and not from a Breit-Wigner fit). Further 
pion- Roper couplings are encoded in and . To analyze the real part 
of the Roper self-energy, one has to calculate a) tree graphs with insertion 
~ c\ , e\ , self-energy diagrams with intermediate b) nucleon and c) Roper 

( 2) 

states and d) tadpoles with vertices from C R . In fact, the graphs of type b) 
require a modification of the regularization scheme due to the appearance 
of the two large mass scales mjv and itir. The solution to this problem - 
assuming m%/m R <C 1 (in nature, this ratio is ~ 1/2.4) - is described in. 14 
As discussed in that paper, the LECs c* and e* can be bounded assuming 
naturalness and by direct comparison with the corresponding pion-nucleon 
couplings: |c*| < 0.5 GeV" 1 , |c^ 3 | < 1.0 GeV" 1 and |ej| < 0.5 GeV- 3 . In 
Fig. 2 an estimated range for the pion mass dependence of the Roper mass 
is presented by taking the extreme values for c* 2 3 and e\, while keeping 
c* x = — 0.5 GeV -1 , gNR — 0.35, gR = 1 fixed. The masses of the baryons 
in the chiral limit are taken to be = 0.885 GeV 16 and rriR — 1.4 GeV, 
respectively. The dash-dotted curve is obtained by setting the couplings 
c* 2 3 , e\ all to zero, and exhibits up to an offset a similar quark mass depen- 
dence as the nucleon result (dotted curve, taken from Ref. 16 ). It should be 
emphasized, however, that the one-loop formula cannot be trusted for pion 
masses much beyond 350 MeV. No sharp decrease of the Roper mass for 
small pion masses is observed for natural values of the couplings. Note that 
the important Arr and Ntttt channels are effectively included through the 
dimension two and four contact interactions, still it would be worthwhile 
to extend these considerations including the delta explicitely. Note further 
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Fig. 2. Quark mass dependence of the Roper mass for different parameter sets c\ = 
—0.5, 3 , ej. The ci are in units of GeV -1 and e\ is given in GeV -3 . and couplings gn = 
1.0, gNR = 0.35. The solid curve corresponds to c| = 1.0, Cg = 1.0, = 0.5, the dashed 
one to = — 1.0, Cg = — 1.0, e? = —0.5 and the dot-dashed one to = = = 0. 
The dotted curve represents the quark mass dependence of the nucleon, see Ref. 16 The 
values of the corresponding LECs are: c\ = — 0.9, C2 = 3.2, C3 = — 3.45,ei = —1.0. 

that the formalism developed in 14 is in general suited to study systems with 
two heavy mass scales in addition to a light mass scale. In this sense, it can 
be applied to other resonances as well, such as the 5n(1535). In this case, 
however, an SU(3) calculation is necessary due to the important r]N decay 
channel. 

4. Intoduction II: Hadronic atoms 

Let us come back to the real world of physical quark masses. In what follows, 
I will discuss a spectacular effect of isospin violation in pionic deuterium, 
which is one particular hadronic atom. More generally, such atoms are made 
of certain hadrons bound by the static Coulomb force. There exist many 
species, e.g. n + ir~ , tt^K^, ir~p, n~d, K~p, or K~d. In these systems, the 
Bohr radii are much larger than any typical scale of strong interactions 
(QCD), so that their effects can be treated as perturbations. These are the 
energy shift AE from the Coulomb value and the decay width T (often 
combined in the complex valued energy shift). Since the average momenta 
in such systems are very small compared to any hadronic scale, hadronic 
atoms give access to scattering at zero energy and thus the pertinent S- 
wave scattering length(s). As it is well known, these scattering lengths are 
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very sensitive to the chiral and isospin symmetry breaking in QCD. In 
fact, hadronic atoms offer may be the most precise method of determining 
these fundamental parameters since theory and experiment can in some 
cases be driven to an accuracy of one or a few percent. On the theoretical 
side, hadronic atoms can be analyzed systematically and consistently in the 
framework of non-relativistic low-energy EFT including virtual photons, see 

17 

e.g. 

5. Isospin violation in pionic deuterium 

Combined measurements of the energy shift and decay width of pio- 
nic hydrogen and the energy shift of pionic deuterium offer an excellent 
test of isospin symmetry and its breaking because these three quantities 
are expressed in terms of two scattering lengths, AE(tt~p) <~ a + + a~ , 
T(ir~p) ~ a~, and ReAE(ir~d) <~ a + + where the ellipses denotes 
three-body effects such as multiple scattering within the deuteron. Here, a + 
and a~ are the isoscalar and the isovector S-wave irN scattering lengths, 
respectively. The Bern group has championed the EFT treatment of pio- 
nic hydrogen, the calculations including strong and electromagnetic isospin 
violation can be found in 18 and 19 for the ground state energy and the 
width, respectively Using this formalism to analyze the data from PSI (as 
reviewed in 20 ) and combining these with the EFT treatment of pionic deu- 
terium in the isospin limit 21 (compare the bands denoted hydrogen energy, 
isospin breaking and hydrogen width, isospin breaking and deuteron, no 
isospin breaking in Fig. 3) one faces a problem - these bands do not in- 
tersect. Since the analysis of pionic deuterium in Rcf. 21 was done in the 
isospin limit, the question naturally arises whether this is the source of 
the trouble? In principle, isospin violation (IV) leads to corrections in the 
bound-state as well as in the pertinent scattering amplitude. From experi- 
ence with pionium, pionic hydrogen and irK atoms (see e.g. Rcf. 22 ), such 
bound-state corrections are expected to be small. On the other hand, al- 
ready in 1977 Weinberg pointed out that IV effects can be unnaturally 
large if the isospin-conserving (IC) contribution is chirally suppressed. 2 
In particular, such an effect is very pronounced in neutral pion scatter- 
ing off nucleons (for an update, see 24 ), but it is very hard to observe. On 
the other hand, the leading order contribution to wd scattering is chirally 
suppressed, Rca^ ~ (a n - p + a w - n ) = 0(p 2 ), despite the fact that a v - p 
and a v - n a re individually of 0(p). Here, p denotes collectively the small 
parameters of CHPT. While this is well-known, nobody has ever system- 
atically investigated IV in ir~d. The leading order IV in pionic deuterium 
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was only analyzed recently in Rcf. 25 To be specific, consider the threshold 
pion-deuteron scattering amplitude: 

Re = Re af d + Aa^ , (4) 

where the IV piece ~ Aa^ appears at the same order as the leading IC 
piece ~ Re a^°J . Also, to this order one has no dependence on the deuteron 
structure. The explicit calculation leads to 

Aa^ = ( 4 ^(l + M /2))- 1 ( ( 5r p + ( 5r„). (5) 

Here, \x = M n + /m p and T Pin are the leading isospin breaking corrections 
to the n~p and ir~n threshold scattering amplitudes. These are given by 

^ ""■y^v ^M-t,^-), ( 6 ) 

where = 92.4 MeV is the pion decay constant, gA = 1-27 denotes the 
axial- vector charge of the nucleon and c\ is a strong and f\ , ji are elec- 
tromagnetic 0(p 2 ) LECs, respectively. Note also that at lowest order in 
CHPT ci is directly related to the value of the pion-nucleon cr-term and 
ji to the proton-neutron mass difference. In the numerical calculations we 
take ci = -0.9^-j GcV~\ 7 f 2 = -(0.97±0.38) GeV -1 . 18 ' 27 Note that the 
errors on the LEC c\ are most conservative. The largest uncertainty in the 
results is introduced by the constant /i , whose value at present is unknown 
and for which the dimensional estimate |/i| < 1.4 GeV -1 has been used. 
Note also, that the hydrogen energy band, which is shown in Fig. 3 corre- 
sponds to the new value of c\ given above. From this, the leading order IV 
contribution to pionic deuterium follows as: 

AaLO = (4 7r (l + v/2))- 1 (5T p + 5T n ) 

with AM 2 = M 2 + — M 2 the squared charged-to-neutral pion mass differ- 
ence. Substituting numerical values for the various low-energy constants, 
which were specified above, one obtains that the correction at 0(p 2 ) is 
extremely large 

A fl LO = _ ( o.oii ™)M- 1 , (8) 

that is Aa^/Rea^ p = 0.42 (central values), using the experimental value 
Re = -(0.0261 ± 0.0005)M" 1 . 26 Moreover, one can immediately see 
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Fig. 3. Determination of the ttN S-wave scattering lengths a + and a~ from the com- 
bined analysis of the experimental data on the pionic hydrogen energy shift and width, 
as well as the pionic deuterium energy shift (details in the text). The cross denoted as 
Beam et al is taken from Ref. 21 The second cross corresponds to the scattering lengths 
given in Eq. (9). 

that the correction moves the deuteron band in Fig. 3 in the right direction: 
the isospin-breaking corrections amount for the bulk of the discrepancy be- 
tween the experimental data on pionic hydrogen and deuterium. Including 
the corrections AaJJ, all bands now have a common intersection area in 
the a + ,a~-plane, see Fig. 3. The resulting values for the irN scattering 
lengths are: 

a+ = (0.0015 ± 0.0022) M" 1 , 

a- = (0.0852 ± 0.0018) M" 1 . (9) 

Further, using the hydrogen energy shift to estimate the LEC /i, we obtain 

h = -2-i^; 2 2 GeV- 1 , (10) 

which is consistent with the dimensional analysis and a recent evaluation 
based on a quark model. 28 Note that the error displayed here does not 
include the uncertainty coming from the higher orders in CHPT and should 
thus be considered preliminary. 
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As we see, the presence of the 0(p 2 ) LECs in the expressions for the 
isospin-breaking corrections leads to a sizeable increase of the uncertainty 
in the output. In order to gain precision, in the fit one might also use those 
particular linear combination(s) of the experimental observables that do 
not contain /i and c\. However, it should be pointed out that such a fit 
imposes much more severe constraints on the data than the uncorrelated 
fit considered above. In fact, applying isospin-breaking corrections only at 
0(p 2 ), we find that the data are still over-constrained in the combined fit. 
From this we finally conclude that to carry out such a combined analysis 
with the required precision, one would have e.g. first to evaluate the isospin- 
breaking corrections with a better accuracy. 

Up to now, we have restricted ourselves to the leading-order isospin- 
breaking correction in CHPT. Calculations at 0(p 3 ) exist only for the 
hydrogen energy shift and yield S e — (—7.2 ± 2.9) • 1CP 218 (using ci = 
(-0.93 ± 0.07) GeV -1 ). The corrections to 0(p 2 ) result are sizable (the 
energy band in Fig. 3 will be shifted further upwards) , but the uncertainty, 
which is almost completely determined by the 0(p 2 ) LECs, remains prac- 
tically the same. On the other hand, consistent studies at 0(p 3 ) imply the 
treatment of the scattering process in the three-body system in the effec- 
tive field theory with virtual photons. To the best of our knowledge, such 
investigations have not been yet carried out, although certain three-body 
contributions at 0(p 3 ) were calculated in the past. 29 It is natural to expect 
that generally three-body terms at 0(p 3 ) should not depend on the addi- 
tional LECs from the two-nucleon sector and hence the extraction of the 
ttN scattering lengths at a high precision is still possible. Of course, these 
arguments can not be a substitute for a rigorous proof in the framework 
of EFT, which in the light of the above discussion, is urgently called for. 
It was also shown in Ref. 25 that the 0(p 4 ) correction which emerges from 
the double-scattering term in the multiple-scattering series is very small, 
Aa doubic scat. = 0.003 Re a^ d p , using the scattering lengths from Eq. (9) 
as input. Of course, such a partial result can only be considered indica- 
tive. Evidently, the systematic analysis of all 0(p 3 ) (and eventually 0(p 4 )) 
corrections should be carried out. 

6. Summary and outlook 

In the first part of this talk, I have considered aspects related to baryon 
CHPT for light quark masses above their physical values. It should be 
stressed that baryon CHPT is a mature field in the up and down quark sec- 
tor and provides unambiguous extrapolation functions for LQCD - certainly 
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more work is needed for the three-flavor case. To carry out the required 
chiral extrapolations, one should keep in mind that different observables 
are linked by general operator structures and the appearing low-energy 
constants (LECs) are universal, which means that they are independent 
of the process considered. Given the present status of LQCD, it appears 
mandatory to perform global fits to observables and constrain the appearing 
LECs by input from phenomenology, whenever available. I have discussed 
two specific examples of the interplay between CHPT and LQCD. A chi- 
ral extrapolation for gA exists now at two-loop accuracy. For pion masses 
< 350 MeV, the theoretical uncertainty related to it is reasonably small. I 
have also provided a chiral extrapolation function for the Roper mass - cer- 
tainly much more work is needed for such excited states from both CHPT 
and LQCD. Evidently, we need more lattice "data" at low quark masses to 
really perform precision studies. 

In the second part of this talk, I returned to the real world (physical 
quark masses) and considered hadronic atoms. These can be systematically 
analyzed in non-relativistic effective field theory including virtual photons. 
We have found a very large isospin-violating effect in pionic deuterium at 
leading order. That there is such an effect is not so surprising because the 
leading isospin-conserving contribution is chirally suppressed. What is sur- 
prising, however, is the actual size of the effect and that it was only found 
recently. Combining this with the information obtained from the analysis 
of the energy shift and width in pionic hydrogen, one is led to a consis- 
tent extraction of the S-wave nN scattering lengths and can furthermore 
determine the electromagnetic LEC f\. Clearly, higher order calculations 
are necessary to reduce the theoretical uncertainty. In this context it is also 
important to stress that it was recently shown that there are only tiny dis- 
persive corrections to Re a^d, see 30 and Hanhart's talk at this conference. 
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